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Abstract 

The flexoelectric (FE) effect provides a linear coupling between 
electric polarization and orientational deformation in liquid crystals. 
It influences many electrooptical phenomena and it is used in some 
bistable nematic devices. A statistical theory of dipole FE polariza- 
tion in biaxial nematic liquid crystals is used to calculate temperature 
dependence of order parameters, elastic constants, and FE coefficients. 
The splitting of the two Meyer FE coefficients and the appearance of 
new FE coefficients is obtained at the uniaxial-biaxial nematic tran- 
sition. The ordering of the splited FE coefficients corresponds to the 
ordering of the splited elastic constants. 

Keywords: biaxial nematic liquid crystals, flexoelectric effect, elastic 
constants. 



1 Introduction 

The biaxial nematic liquid crystals were first predicted by Preiser [1], who 
showed that molecules with shapes that deviate from cylindrical symmetry 
could possess a nematic phase with three distinct optical axes. A biaxial 
nematic phase was first observed in a lyotropic mixture by Yu and Saupe [2] 
in 1980 but the existence of a thermotropic biaxial system was not certain for 
many years. Several reports appeared in 2004 on thermotropic liquid crystals 
formed by side-chain polymers [3j, bent-core molecules |4l[5], and tetrapodes 
[6]. The number of new biaxial systems is constantly growing. 



1 



Many theoretical papers [7] and computer simulations [S] show that the 
molecular shape and pair interaction biaxiality are important for the biaxial 
phase to exist. However, very often real systems favour packing in the smec- 
tic or crystalline biaxial phases. It is a challenge for the theory to find factors 
responsible for absolute stability of the biaxial nematic phase. It was shown 
that fluctuations in molecular shape can influence the biaxial nematic phase 
stability jl]. The motivation for this search ranges from purely academic 
interest to the potential usage of these materials in faster displays, where 
in principle the commutation of the secondary director should give lower re- 
sponse times compared to the conventional twisted nematic and ferroelectric 
smectic devices. 

A static electric field imposed on a nematic liquid crystal have many phys- 
ical effects, but the most important are two of them. One is connected with 
the anisotropy of the dielectric constant. The second effect is the appearance 
of the spontaneous polarization in a deformed liquid crystal; this is called 
the fiexoelectric effect. Conversely, an electric field may induce distortions 
in the bulk. In 1969 Meyer showed that it is a steric effect due to the shape 
asymmetry of polar molecules fTU]. In case of nonpolar molecules the FE 
effect originates from a gradient of quadrupole moment density [llj. 

The two FE coefficients were introduced by Meyer for splay and bend 
distortions of the uniaxial nematic phase [lOj. Recently, a statistical theory 
for the dipole FE polarization was derived in the case of the biaxial nematic 
phase composed of the molecules [12j. There are six splay-bend defor- 
mations of the biaxial nematic phase and thus six FE coefficients are defined, 
but only five of them are independent. General microscopic expressions for 
the FE coefficients involve the one-particle distribution function and the po- 
tential energy of two-body short-ranged interactions. 

The FE effect has a large infiuence on many phenomena in liquid crys- 
tals: electrooptical phenomena and defect formation, for instance. It plays 
a key role in some device applications. Fiexoelectric switching is important 
in bistable displays [13]- [15]. Fiexoelectric coupling in chiral and twisted ne- 
matic crystals [TB] leads to a linear rotation of the optic axis and also leads 
to device applications [T7]. Fiexoelectric coupling in smectic liquid crystals 
has been shown to stabilize helical structures [18j . 

The purpose of this study is to calculate the values of the FE coefficients 
and other material parameters for model systems. The proper form of the 
interaction potential energy allows us to calculate the temperature depen- 
dence of the order parameters, the FE coefficients, and the elastic constants. 



2 



The uniaxial and biaxial nematic phases are considered. 

2 Description of the system 

Let us consider a set of molecules contained in a volume V, at the tem- 
perature T. The molecules are rigid blocks {C2V symmetry) with three trans- 
lational and three rotational degrees of freedom. It is assumed that the 
molecules interact via two-body short-range forces that depend on the dis- 
tance between the molecules {u = f2 — fi = uA) and their orientations de- 
scribed by the three Euler angles R = (0, 6, ip) or three orthonormal vectors 
(/, m, n). 

The microscopic free energy F for the system is given by 

F = Fent + Fint, (l) 



/3Fent = / drdRG{f, R){ln[G{f, R)A] - 1}, (2) 
PF,nt = -^J dndRidf2dR2G{n,R,)Gif2,R2)fu. (3) 

Here /12 = exp(— /3$i2) — 1 is the Mayer function, $12 the potential en- 
ergy of interactions, dR = d(f)d6 sin 6 dip, f3 = l/lksT), and A is related to 
the ideal gas properties. The one-particle distribution function G has the 
normalization 

J dfdRG{f, R) = N. (4) 

The equilibrium distribution G minimizing the free energy ([T]) satisfies the 
equation 

ln[G'(fi, i?i)A] - J df2dR2G{r2, i?2)/i2 = const. (5) 

For the homogeneous phase G = Gq does not depend on the position and it 
has the form [T9j 

Go{R) = Go{l- L,l- N,n ■ L,n- N), (6) 

were the orthonormal vectors (L, M, A^) define the biaxial nematic phase 
axes. In practice we characterize the alignment not through the full function 
G, but by some numerical parameters - order parameters. In the case of 
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the biaxial nematic phase the main order parameters are the orientational 
distribution averages of the following four functions [20j: Fqq , Fq2\ , and 
f!^ . We note that there are other notations [22j. In the uniaxial nematic 
phase the functions F^ and F^^ have nonzero averages only. The molecule 
alignment can be also described by order tensors Q which are often calculated 
for computer simulations [23] where the molecular and laboratory axes must 
be distinguished. The order tensors are defined as 

Q% = i3Lh-6.^)/2, (7) 
Q'^P = {3m^mp-6ap)/2, (8) 
Q^^ = i3n^n^-6^p)/2. (9) 



3 Elastic deformations of the phase 

Orientational ordering of biaxial nematics is usually described by the thee 
orthonormal vectors 

L = RiaCa, M = R2aea, N = RsaGa- (10) 

In the homogeneous phase the vectors (L, M, A^) are constant in space, but 
in a deformed phase they depend on the position in space. In a continuum 
approach the distortion free-energy density is obtained as an expansion 
about an undistorted reference state with respect to gradients of the vectors 
{L, M, N). The form of the fd can be derived in many alternative ways but 
we use the form presented by Stallinga and Vertogen [21j (the surface terms 
are neglected) 



^K^niiDii)' + ^K^2i2{Du? + ^^^1313(^13)' 
X 2 1 2 1 2 

+ 2-^2121 (-D21) + 2-^2222(-D22) + 2-^2323(-D23) 

+ 2-^3131 (-Dsi)^ + --^"3232 (.032)^ + 2-^3333(-D33)^ 
+-^1122-Dii-D22 + -^1133-Dii-D33 + K2233D22DS3 
+_K'l221-Dl2-D21 + -^"1331-013-031 + -^^2332 -023-032- (H) 

Dij = -^jklRiaRkl3daRll3 ■ (12) 
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Microscopic expressions for the elastic constants Kiji^i were derived in [12] 
and it was shown that there are 12 independent bulk constants because 

K122I = -^1122, Ki32i = K1133, i^2332 = -^^2233- (13) 

4 Flexoelectric polarization 

Liquid crystalline phases often consist of polar molecules but in homoge- 
neous nematic phases the average polarization is zero. On the other hand, a 
phase distortion can produce a polarization and this is called the FE effect. 
In a continuum approach the FE polarization of the biaxial nematic phase 
depends on the spatial derivatives of the vectors (L, M, A^) [T5] 

Pa = '^{SiiRiadfiRili + huRijidjiRia). (14) 
i 

The parameters Su and ha, {i = 1, 2, 3) are not unique because if we add 
any constant to all of them, the polarization will not change. The physical 
FE coefficients (i = 4, . . . , 9) are 

'^4 = S33 — 611, 05 = S22 — bii, = S33 — 622, 

aj = sn - 622, as = S22 - ^33, ctg = su - 633. (15) 

The coefficients satisfy the identity 

04 — 05 — ae + ay + ag — ctg = 0. (16) 

Deformations of the biaxial nematic phase connected with the FE effect are 
given in Tabled] In the case of the uniaxial nematic phase the FE polarization 
has the form 

P„ = eiN^dpNp + e^NpdpN^. (17) 
Let us define the molecule electric dipole moment as 

A^Q = /iiL + Ai2"^Q + yWs^a- (18) 

In the case of the molecular iteractions described below, the FE coefficients 
can be expressed as follows 

04 = dudRidR2fi2Go{Ri)fi3ni,{-u^){U2z-W2x), (19) 
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05 = 


1 dudRidR2fi2Go{Ri)fi3niy{-Ux)U2y, 


(20) 


oe = 


J dudRidR2fi2Go{Ri)iJ,3ni^UyW2y, 


(21) 


ay = 


J dudRidR2fi2Go{Ri)fi3ni^UyU2y, 


(22) 


as = 


J dudRidR2fi2Go{Ri)iJ,3niy{-u^)W2y, 


(23) 


ag = 


J dudRidR2fi2Go{Ri)n3ni^u^{U2z -W2x), 


(24) 



where it is assumed that n defines the molecule G2 axis and 

Ua = diGoL + d^GoUa, Wa = d2Gola + d^GoUa. (25) 

5 Results 

We performed our calculations for the square-well potential energy of the 
form 

{+00 for (u/a) < 1, 
-e for 1 < (u/a) < 2, (26) 
for (u/a) > 2, 

were e is the depth of the well and a depends on the molecule orientations 
and on the vector A 

a = ao + (Ti(A-ni-A-n2) + a2 [(A-fii)2 + (A-n2; 



+ 0-3 



(A ■ + (A • I2Y (27) 



The parameter ctq defines the length scale, cxi defines the FE term, (T2 and 
(T3 define biaxial nematic terms. We used the density NVmoi/V = 0.1, the 
molecule volume Vmoi was estimated from the mutually excluded volume. 
The FE coefficients were expressed in Hi/ctq, the elastic constants in e/cxo, and 
the temperature in e/A;^. The parameters cTj are given in Table [2l The two 
physical systems are considered that consist of wedge-shaped and banana- 
shaped molecules. 

5.1 Wedge-shaped biaxial molecules 

In the system of wedge-shaped biaxial molecules, on decreasing the temper- 
ature we meet the first order transition from the isotropic to the uniaxial 
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nematic phase at T/^r = 0.618 and the second order transition to the biaxial 
nematic phase at T/vb = 0.401. The temperature dependence of the order 
tensors is presented in Fig. [H The temperature dependence of the elastic 
constants and the FE coefficients are presented in Figs. [2] and [3l respectively. 

The values of Q"" show that long molecule axes are directed along the 
Z axis in the whole nematic region, whereas the values of Q^'^. reveal the 
alignment of short molecule axes along the X axis and it is enhanced in 
the biaxial nematic phase. The splay elastic constant Ki splits into -ft"i2i2 
and -ftr2i2i- The bend elastic constant splits into -ft'3232 and -ft'sisi. Note 
that the equality Ki = is accidental and results from neglecting order 
parameters F^{) with j > 2. The splay FE coefficient Ci splits into 04 and 
(a4 > Og > 0). The bend FE coefficient —63 splits into and Og (0 > ag > 
ag). The coefficients 05 and ay are small and almost always negative. 

5.2 Banana- shaped biaxial molecules 

In the system of banana-shaped biaxial molecules, on decreasing the tem- 
perature we meet the first order transition from the isotropic to the uniaxial 
nematic phase at Tin = 0.595 and the second order transition to the biaxial 
nematic phase at Tnb = 0.382. The temperature dependence of the order 
tensors is presented in Fig. HI The temperature dependence of the elastic 
constants and the FE coefficients are presented in Figs. [S] and El respectively. 

According to the values of Q'j^ long molecule axes are directed along the 
Z axis in the whole nematic region, whereas the values of show the 
alignment of short molecule axes along the X axis and it is enhanced in the 
biaxial nematic phase. The behaviour of the elastic constants is similar to 
the case of the wedge-shaped molecules because the FE term is small in both 
cases. The bend FE coefficient —63 splits into ag and ag (0 > ag > ag). The 
splay FE coefficient ei is smaller then 63 and it splits into a4 and a^. The 
coefficients as and ay are again small but comparable with a4 and Qq. The 
sign of some coefficient can change on changing the temperature. 

6 Conclusions 

In this paper, the statistical theory was used to study the temperature de- 
pendence of the order parameters, elastic constants, and FE coefficients of 
biaxial nematic liquid crystals. In order to calculate these macroscopic pa- 
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rameters one needs the one-particle distribution function and the potential 
energy of molecular interactions. The two physical systems were considered. 
The splittings of the FE coefficients and the elastic constants were obtained at 
the uniaxial-biaxial nematic transition. New small FE coefficients appeared 
at the transition. The ordering of the splited FE coefficients corresponds to 
the ordering of the splited elastic constants. 

The FE coefficients were proportional to the dipole moment component 
parallel to the molecule C2V symmetry axis. This was the result of the inter- 
actions potential symmetry. The beaviour of the main FE coefficients, ei for 
the wegde-shaped molecules and 63 for the banana-shaped molecules, is clear 
and it is in the agreement with previous studies On the other hand, 

it seems that other FE coefficients should be interpreted with caution. It is 
possible that higher order parameters can have a significant contribution. 

At present stage, the direct comparison between the theory and the ex- 
periment in not possible because to our knowledge the FE coefficients have 
not been measured for the biaxial nematic phase. What is more, even for 
the uniaxial nematic phase the experimental data are still scarce and some- 
times contradictory [25] . However, when biaxial nematic phases become more 
widespread, the presented theory will be helpful in practical applications. 
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Table 1: Deformations of the biaxial nematic phase connected with the FE 
effect. The corresponding elastic constants and the FE coefficients are given, 
the values for the uniaxial nematic phase are in parentheses. 



Deformation Elastic constant FE coefficient 



— =? =? 

N splay, L bend -?^i2i2 (-^i) 04 (ei) 

M splay, L bend -^1313 (0) 05 (0) 

N splay, M bend ^^2121 (-^1) ^6 (ci) 

L splay, M bend K2323 (0) ay (0) 

— * — * 

M splay, N bend -f^sisi (Ks) as (—63) 

L spltw. bend A'3232 ij'^'s) ^'9 (^'':-!) 



Table 2: Parameters CTj used in calculations. 



Molecules 


(Ji/(Jo 






Long axis 


Short axis 


C2 axis 


wedgc-likc 


0.2 


0.5 


-0.4 


n 


I 


n 


banana-like 


0.2 


-0.4 


0.5 


— * 


n 


n 
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Figure 1: The temperature dependence of the order tensors for wedge-shaped 
biaxial molecules. 
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T 

Figure 2: The temperature dependence of the elastic constants for wedge- 
shaped biaxial molecules. The squares, triangles, and circles denote defor- 

— * — * — * 

mations with N splay, bend, and constant, respectively. The empty 
(filled) symbols indicate the larger (smaller) parameter. 
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Figure 4: The temperature dependence of the order tensors for banana- 
shaped biaxial molecules. 
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